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Fine structures in the tunneling spectra of superconductors have been widely used to identify 
fingerprints of the interaction responsible for Cooper pairing. Here we show that for scanning 
tunneling microscopy (STM) of Pb, the inclusion of inelastic tunneling processes is essential for the 
proper interpretation of these fine structures. For STM the usual McMillan inversion algorithm of 
tunneling spectra must therefore be modified to include inelastic tunneling events, an insight that is 
crucial for the identification of the pairing glue in conventional and unconventional superconductors 
alike. 

PACS numbers: 74.55.+V, 74.81.Bd, 74.25.Jb, 74.25.Kc 


Conventional superconductivity is caused by the at¬ 
tractive interaction between electrons near the Fermi en¬ 
ergy mediated by phonons [I]. This leads to the for¬ 
mation of a gap 2A in the single particle density of 
states (DOS) of the electrons, and to quasi-particle peaks 
above and below the gap m n. Eliashberg extended 
the BCS theory to the limit of larger dimensionless 
electron-phonon coupling constants A, included a real¬ 
istic electron-phonon coupling and the detailed structure 
of the phonon spectrum [3]. As a consequence, the quasi¬ 
particle peaks near the Fermi surface are modified due to 
the interaction with phonons, leading to fine structures in 
the electronic DOS near the peaks of the Eliashberg func¬ 
tion a 2 F(co) shifted by A. a 2 F(to) is the phonon DOS 
F(co), weighted by the energy dependent electron-phonon 
coupling strength a 2 (oj). These fine structures are due 
to the excitation of virtual phonons (see Fig. 1). Ex¬ 
perimentally, these fine structures in the electronic DOS 
have been detected with electron tunneling spectroscopy 
on planar junctions rang. In the pioneering work of 
McMillan and Rowell dl. the Eliashberg function could 
be reconstructed from the superconducting DOS by an 
inversion algorithm taking into account the interaction 
of electrons and virtual phonons. This method has been 
used to identify fingerprints of the phononic pairing glue 
in the electronic spectrum and thus to determine the pair¬ 
ing mechanism leading to superconductivity nn nu. it 
counts as a hallmark of condensed matter physics. 

An alternative way to determine the Eliashberg func¬ 
tion is to measure the energy dependence of the scattering 
of electrons with real phonons in the normal state using 
inelastic tunneling spectroscopy (ITS) [TOUT] , see Fig.l. 
This method is more direct, as the second derivative of 
the tunneling current I with respect to the bias voltage 
U is, under rather general assumptions, directly propor¬ 
tional to a 2 F(oj) [18]. Recently, this method has been 
combined with scanning tunneling microscopy (STM) to 
obtain local information on the Eliashberg function of Pb 
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FIG. 1: Illustration of the inelastic tunneling processes from 
a sharp tip (right) into a superconductor (left) in real space. 
Filled states are shown in red colour with energy along the 
vertical axis. The inelastic tunneling process is accompanied 
by the excitation of real phonons (green). 


on a Cu(lll) substrate [T5] . 

In this work, we determine experimentally and ana¬ 
lyze theoretically the tunneling conductance of Pb that 
is affected by the coupling to real phonons via inelas¬ 
tic tunneling and virtual phonons via many-body renor¬ 
malizations. Comparing the two approaches to deter¬ 
mine a 2 F(ui ) on the same sample with the same tip of 
a low temperature STM, we show that interpreting tun¬ 
neling spectra of superconductors via the McMillan in¬ 
version algorithm (and thus solely by its elastic contribu¬ 
tion) can be an incomplete description. We demonstrate 
that inelastic contributions to the tunneling current can, 
in general, be of the same order as the elastic contri¬ 
bution. We show that we can understand experimental 
STM data from Pb tunneling in the normal and super¬ 
conducting state, taking into account both elastic and 
inelastic tunneling processes. The combined analysis of 
elastic and inelastic tunneling processes is important to 
correctly identify fingerprints of the relevant interactions 
in the electronic DOS and to identify the pairing glue 
for superconductivity. This is essential for conventional 
superconductors, such as Pb, but is expected to be even 
more important for unconventional pairing states, where 
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an electronic pairing interaction is expected to funda¬ 
mentally change its character upon entering the super¬ 
conducting state. 

We start with experimental data for STM measure¬ 
ments on lead. Measurements were performed with a 
home-build Joule-Thomson low-temperature STM (JT- 
STM) [28] at temperatures about 0.8 K. The JT-STM 
contains a magnet which allows to suppress supercon¬ 
ductivity. In order to ensure that there is no significant 
inelastic signal of the tip at \U\ < 15 mV, we use a chemi¬ 
cally etched tungsten tip, known to have a weak electron- 
phonon coupling [29]. The highly n-doped Si(lll) crys¬ 
tals were carefully degassed at 700 °C for several hours 
and then flashed to 1150 °C for 30 seconds to remove 
the native oxide. Lead was evaporated at room temper¬ 
ature from a Knudsen cell with a nominal thickness of 
19 monolayers (ML). After deposition the samples were 
immediately transferred to the cryogenic STM. In agree¬ 
ment with previous studies flat-top, wedge-like 

islands of local thickness around 30ML were observed 
(see Fig. [2]), i.e. extended 3D islands appear on top 
of a metallic wetting layer (WL) [o5| . The islands are 
Pb single crystals with their (111) axis perpendicular to 
the substrate [2D 133 HO- The first (second) derivative 
of the tunneling current dl/dU ( d 2 I/dU 2 ) of the islands 
was measured using a lock-in amplifier with a modulation 
voltage of ITmod = 439 /rV. 



FIG. 2: STM topography of Pb on Si(lll) (300 x 
175 nm 2 , IV, InA). The thickness of the island was de¬ 
termined to x ~ 30 monolayers. 



FIG. 3: Second derivative d 2 I/dU 2 ~ a 2 Ftun(eiy) measured 
in the normal conducting state (T = 0.8 K, B = 1 T). 


clearly revealing the two characteristic phonon peaks that 
are also seen in the Eliashberg function a 2 F(co) deter¬ 
mined by Ref. m Below we show explicitly that in 
the normal state d 2 I/dU 2 is proportional to a 2 F tun (e U). 
These peaks at U = 4.05 mV« uj t and U « 8.3 mV« oji 
(FWHM 7 1 = 1.076meV and 7 1 = O.GOmeV) coincide 
with the energies of the transversal and longitudinal van 
Hove singularities in the phonon DOS of lead [ 23133 . 
The additional peak at U « 12.5 mV can be explained 
by tunneling processes via two-phonon emission [56] . 

The key implication from Fig. [3] for the superconduct¬ 
ing state is, however, that we must include inelastic con¬ 
tributions to the superconducting tunneling spectrum in 
a consistent fashion. Before we present our experimental 
data of the superconducting state, we summarize the the¬ 
oretical description of the tunneling conductance in the 
superconducting state including inelastic contributions. 

The Hamiltonian H = Ho + Ht used in our analysis of 
the combined substrate and tip consists of free electrons 
in the tip and electrons interacting with phonons in the 
substrate (we set H= 1 ): 


n 0 


e p C p.a C P,cr + e k C \,cr C k,cr + W q,P a \,p. a q,P 

p,cr k,c t q,p 

I 7vP ^ ' a k-k ',/Afc — fc',M ’ (1) 

v Ks 


While the electrons in the ss 30ML Pb film on Si 
have quantized k z leading to the flat island growth, the 
phonon DOS of the finite thickness films is rather similar 
to that of bulk Pb as indicated by first principle calcula¬ 
tions [35t 136] . As a first measurement, we therefore deter¬ 
mine a 2 F tun (uj) of lead directly with ITS in the normal 
state. Pb islands were forced to the normal state by ap¬ 
plying a magnetic field of IT normal to the sample plane. 
Since the sample is in the normal state, no renormaliza¬ 
tion of the BCS density of states near the Fermi energy 
due to virtual phonons arises. Thus renormalization ef¬ 
fects by virtual phonons can be neglected in d 2 I/dU 2 and 
experimental features in d 2 I/dU 2 correspond to inelastic 
tunneling. Fig.[3]shows the measured d 2 I/dU 2 spectrum 


Here, 4> q fl = a 9iM + a*_ q /( is proportional to the lat¬ 
tice displacement, where a q is the the phonon anni¬ 
hilation operator for momentum q and phonon-branch 
and with dispersion oj q c' k / p a are the electron an¬ 
nihilation operators for the two subsystems: The tip 
(quasi-momentum p, dispersion ej and volume Vr) and 
the superconductor (quasi-momentum fe, dispersion ej? 
and volume Vs). For the latter we include the electron- 
phonon coupling af.-k',p that gives rise to superconduc¬ 
tivity. The electron-phonon interaction in the tip is as¬ 
sumed to be small. In addition, the tunneling part of 
the Hamiltonian includes elastic and inelastic tunneling 
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processes [HI1S]: 

1 


y/VrVs J 2 T k.P C {,a C P,° + h - C ' > ( 2 ) 

Tk,p = T k p + -^= r ^k,p,q,p a q,P ( t > q,p + • 


9,M 


The first term of the tunneling amplitude 7), p describes 
the elastic tunneling part, the second term corresponds 
to electron transitions via the emission/absorption of 
phonons, see Fig. [Tj It is proportional to the bulk 
electron-phonon coupling a q>tl [18]. There can also be 
processes with a higher number of phonons, which will 
be discussed later. 

In order to determine the tunneling current we assume 
that the DOS of the tip is constant ~ and that 

the tunneling amplitudes are independent of momenta 
and phonon branches Tjj = t e and Tj, = t\ which 
is a reasonable approximation for STM [2l]. Then, to 
leading order in t e , the differential conductance gives the 
well known result 


e/m dP 

" {u)= w = 


dui n' F (u) + eU)Ds(co). (3) 


In the limit that T is smaller than the electronic en¬ 
ergy scales, the conductance is just proportional to the 
the normalized electron DOS vs(uj) = I'g (w)/z7y, where 
i/g is the normal state DOS of the superconductor at 
the Fermi level. The conductance constant is given by 
<To = Aite 2 \t e \ 2 and n F is the Fermi function. In 
the normal state, z/g( u) is essentially constant for small 
applied voltages and the second derivative of the elastic 
current vanishes, as discussed above. In the supercon¬ 
ducting state, the opening of the superconducting gap 
and the excitation of virtual phonons lead to the men¬ 
tioned fingerprints of superconductivity and the pairing 
glue in the elastic tunneling spectrum. Below we deter¬ 
mine these structures from the solution of the nonlinear 
Eliashberg equations for given a 2 F(to) and compare with 
our STM experiments. 

The inelastic contribution to the differential conduc¬ 
tance cr 1 (C7) = c jjj due to the excitation of single real 
phonons is for U > 0 given by the convolution 

Ui p poo 

a'(U)=cr o U Q / dua 2 F^eU + u)v s (u))n F (u). 

\t e \ Vg < 7—00 

(4) 


in the limit that the thermal phonons can be ne¬ 
glected T <C ojd- The function a 2 F^ n (x) = 
~ f-oodyaPFtunivWpiy - x) is a thermally 
broadened version of the weighted phonon DOS 

a 2 F tun (uj) = y§ E g>/i \ a q F \ 2 “ u q,n) that is 

closely related to the Eliashberg function a 2 F(co) = 


pp/j E k,k'.p \ a k-k’,p\ 2 - u) k - k ' tli )6{e%)6(e%,). Both 

have similar features but can differ in fine-structure and 
amplitude. 

The result © is the generalization of the current in 
the normal state, where 3 ^-| NC ~ sign(Z7) a 2 F tun (e |C/|) 
is proportional to the weighted DOS of the phonons (or 
other collective excitations of the system), see Ref. [TS] 
EDI EE, M- If naturally explains the results of Fig. [5] 
or the recent STM measurements on Pb [ID] . Our mea¬ 
surement further allows for an estimate of the inelastic 
tunneling amplitude t 1 ss t e /D , which is inversely pro¬ 
portional to the characteristic energy scale of the off- 
shell electrons involved in the tunneling process. The 
normal state elastic conductance cr e (U) ss op is not en¬ 
ergy dependent for the applied biases U and we empha¬ 
size that all spectra within this paper are normalized to 
(Jo = <t(0 ) = cr e (0) to point out the existence of inelas¬ 
tic tunneling contributions. The change in the conduc¬ 
tance from 0 to 10 mV seen in Fig. EH is purely due 
to the inelastic tunneling. This leads to the condition 
CT^lOmV) « 12% (T 0 , where er(0) = op is the purely elas¬ 
tic contribution at zero bias. Using the widely accepted 
Eliashberg function a 2 F(oj) and the experimental DOS 
for lead E3, we can estimate for the characteristic off- 
shell electronic energy to be D ss 240 meV. Below, we 
will see that elastic and inelastic contributions to the fine- 
structure turn out to be comparable in magnitude. 

In the superconducting state, the inelastic contribu¬ 
tion Eq. © has its major contribution slightly below the 
energy of the phonon peaks shifted by the gap A. Since 
inelastic tunneling opens additional channels to the con¬ 
ductance, it will lead to positive contributions to d 2 I/dU 2 
at positive bias. Elastic contributions are of opposite sign 
(see ©)• Thus, pronounced peaks in the second deriva¬ 
tive of the tunneling current due to real phonons are fol¬ 
lowed by dips of same amplitude due to virtual phonons 
(for details see discussion of a single phonon mode in the 
Supplementary Material). As we will see below, we find 
exactly these features in the tunneling current for the 
STM experiment on lead. 

Tunneling processes with a higher number of excited 
phonons will give similar terms as in © with higher con¬ 
volutions of the Eliashberg-function such as a 4 F 2 un (w) = 
f duj' a 2 F tun (u) — uj')ar F tun (u)') and one can formally ab¬ 
sorb this contribution in a redefinition of a 2 F tun (see Sup¬ 
plementary Material). 

Without magnetic field, the islands are in the super¬ 
conducting state. As the local thickness of the islands 
(30ML « 10 nm) is significantly smaller than the bulk 
coherence length of lead (83 nm [35]), the superconduct¬ 
ing gap is not fully developed 30.1311 139H43] , which im¬ 
plies that the spectral weight of the coherence peaks is 
accordingly smaller (see Fig. EH ). Besides the Bogoli- 
ubov quasiparticle peak one clearly observes fine struc¬ 
tures in the spectrum of the conductance around U ~ 
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FIG. 4: Differential conductance dl/dU in the normal (a) and 
superconducting (b) state measured on the island marked by 
arrow in Fig. [2] The curves are normalized to the zero bias 
conductance <r(0) in the normal state. 


5.3 mV and U ~ 9.4 mV. These energies correspond to 
the van Hove singularities in the phonon DOS F(u>) of 
lead shifted by the gap A ss 1.2 meV clearly indicat¬ 
ing electron-phonon interaction induced effects. Further¬ 
more, the typical w/vG ; 2 — A 2 behavior in the BCS DOS 
is altered by the emergence of inelastic contributions at 
biases Vq > 5 mV . This is in contrast to previous mea¬ 
surements on planar tunneling junctions of lead [MS], 
where these inelastic contributions were about one order 
of magnitude smaller m than in our present experiment. 
The reason is that inelastic tunneling events are enhanced 
in STM geometries, if compared to planar tunneling junc¬ 
tions, because the momentum conservation for momenta 
parrallel to the surface are less restrictive [Si- 
Let us now investigate the second derivative of the tun¬ 
neling current in the superconducting state, which is sig¬ 
nificantly more sensitive to the fine structure induced by 
the electron-phonon interaction. For the theoretical spec¬ 
trum, we first use a parameterization of the a 2 F (w), fi* 
from McMillan and Rowell nuns to solve the Eliash- 
berg equations numerically [45] to obtain the lead DOS 
i's(w) in the superconducting state. The elastic contri¬ 
bution to the second derivative is then easily calculated 
using Eq. ©• For the inelastic contribution we use the 
a 2 Ft un (w) function (without the negative dip at small 
voltages U < 2 mV that comes from a zero bias anomaly) 
and the calculated DOS J's(w) to determine the con¬ 
volution in Eq. ([ 6 ]), where the usage of the measured 
& 2 Ft un (aj) function automatically includes two-phonon 
processes and yields the correct amplitude for the inelas¬ 
tic tunneling current. Note, that the rapid wiggels on top 
of the calculated inelastic curve are due to noise of the 
input data of the calculations, i.e. the experimental in¬ 
elastic spectrum in the normal state. This noise is caused 
by residual mechanical vibrations on the level of 300 fm. 
When convoluting the noisy experimental spectra with 
the DOS for the calculation of the inelastic contribution 
in the superconducting state, certain frequencies of the 
noise are amplified and show up as small wiggels. Finally, 
we convoluted the elastic part m with a Gaussian dis¬ 
tribution (standard deviation a = 310/xeV corresponding 


to an energy resolution of 744/reV), describing the exper¬ 
imental broadening due to the modulation voltage of the 

lock-in detection [ra¬ 
in Fig. [5] we compare the experimental data with the 
theoretical prediction of the elastic, inelastic and total 
contributions of the second derivative of the current. 
The experimental data show peak-dip features around 
the zero axis at positions that correspond to the char¬ 
acteristic longitudinal and transversal phonon energies 
w t /i shifted by the gap A w 1.2 meV. For both fea¬ 
tures there is a positive peak at ~ co t ,i + A — 7 ^ of the 
same magnitude as the corresponding dip at u>t,i + A, 
where 7 tj ; are the half-widths of the phonon peaks ob¬ 
served in Fig. [3] This is in contrast to the theoretical 
elastic d 2 I e /dU 2 ~ u'(—eU) curve, which only shows the 
typical dips around A + w t /i predicted by the Eliash- 
berg theory. We note that conventional Eliashberg the¬ 
ory can also have positive peaks, but the following dip 
will always be significantly more pronounced (see also 
Fig. 4 in the Supplementary material). Therefore, the 
observed peak-dip features cannot be explained by pure 
elastic tunneling. However, the measured spectrum both 
in the normal and in the superconducting state can nat¬ 
urally be explained when we combine inelastic and elas¬ 
tic contributions. As can be seen, the total theoretical 
conductance d 2 I tot /dU 2 consisting of elastic and inelas¬ 
tic channels fits the experimental peak-dip features much 
better at the correct energies. 

In summary, we demonstrated experimentally and the¬ 
oretically that in normal conducting Pb islands it is pos¬ 
sible to directly measure the collective bosonic excita¬ 
tion spectrum, here phonons, using STM. In the normal 
conducting state, the obtained d 2 I/dU 2 spectra is pro¬ 
portional to the weighted phonon DOS a 2 F tun (uj) and 
higher convolutions thereof. This is different in the su¬ 
perconducting state of Pb. Here, the obtained second 
derivative d 2 I/dU 2 = d 2 I e /dU 2 + d 2 I'/dU 2 spectra are 
a composition of elastic and inelastic tunneling processes 
with fine structures in the same energy regime. While 
the elastic part shows phonon features coming from self 
energy corrections (exchange of virtual phonons) that 
appear mainly as dips in the second derivative of the 
tunneling current, the inelastic part shows features of 
a 2 i 7 ’ tun (w) shifted by the superconducting gap A giving 
rise to additional peak features of the same amplitude 
at lower energies (excitation of real phonons). A rather 
unique signature of these inelastic contributions are peak- 
dip features in d 2 I / dU 2 around zero at A + w p h in the 
superconducting state. Those cannot be explained by 
only taking into account the elastic part d 2 I e /dU 2 . For 
this reason, the neglect of inelastic processes in STM ex¬ 
periments in general not justified. Hence, when analyz¬ 
ing STM tunneling spectra via the McMillan inversion 
algorithm mm, that gives the purely elastic contribu¬ 
tion, one should carefully subtract the inelastic contri¬ 
butions from the experimental tunneling current. Other- 
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FIG. 5: Comparison of experimental data (red) and theoretical prediction in the superconducting state: Calculated elastic 
(blue), inelastic (green) and total (black) contribution to d 2 I/dU 2 (the elastic current is convoluted with a Gaussian function 
with standard deviation a = 0.31 meV simulating the experimental broadening due to the modulation voltage of the lock-in 
technique). Characteristic peak-dip features around the zero axis can only be explained taking into account elastic and inelastic 
channels (7 tot = 7 el + 7 inel ). 


wise grossly incorrect conclusions about the pairing glue 
would be deduced from the tunneling spectrum. 

Having found out experimentally and theoretically 
how elastic and inelastic tunneling can be disentangled 
for STM in conventional superconductors, the approach 
can be generalized to the investigation of corresponding 
bosonic structures in high temperature superconductors 
such as cuprates and iron pnictides in the future. A cru¬ 
cial difference to the phononic pairing glue is that in 
case of electronic pairing, the bosonic spectrum under¬ 
goes dramatic reorganization below T c in form of a sharp 
resonance in the dynamic spin excitation spectrum [iBh 
EB. Our results imply that great care must be taken 
in the proper interpretation of the tunneling spectra of 
these systems and that real and virtual bosonic excita¬ 
tions must be disentangled in a fashion similar to our 
analysis for lead. 

The authors acknowledge funding by the DFG under 
the grant SCHM 1031/7-1 and WU 349/12-1. 


Ill J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Phys. Rev. 
108 , 1175 (1957). 

[2] I. Giaever, Phys. Rev. Lett. 5, 147 (1960). 

[3] J. Nicol, S. Shapiro, and P. H. Smith, Phys. Rev. Lett. 
5, 461 (1960). 

[4] G. M. Eliashberg, Sov. Phys. JETP 11 , 696 (1960). 

[5] I. Giaever, H. R. Hart, and K. Mergele, Phys. Rev. 126 , 
941 (1962). 

[6] J. R. Schrieffer, D. J. Scalapino, and J. W. Wilkins, Phys. 
Rev. Lett. 10 , 336 (1963). 

[7] J. M. Rowell, A. G. Chynoweth, and J. C. Phillips, Phys. 
Rev. Lett. 9 , 59 (1962). 

[8] J. M. Rowell, P.W. Anderson, and D.E. Thomas, Phys. 
Rev. Lett. 10 , 334 (1963). 

[9] W. McMillan and J. Rowell, Superconductivity Vol. 1, 
ed. R. D. Parks (Dekker, New York, 1969) pp. 561-611. 


[10] I. Giaever, Science 183 , 1253 (1974). 

11 W. L. McMillan and J. M. Rowell, Phys. Rev. Lett. 14, 
108 (1965). 

[12] D.J. Scalapino, J. W. Wilins, Phys. Rev. 148, 263-279 
(1966). 

[13] J.P. Carbotte, Rev. Mod. Phys. 62, 1027 (1990). 

[14] A. Leger and J. Klein, Phys. Lett. 28A, 751 (1969).. 

[15] J. M. Rowell, W. L. McMillan, and W. Feldmann, Phys. 
Rev. 180 , 658 (1969). 

[16] W. Wattamanuik, H. Kreuzer, and J. Adler, Physics Let¬ 
ters A 37 , 7 (1971). 

[17] J. Klein, A. Leger, M. Belin, D. Defourneau, and M. J. 
L. Sangster, Phys. Rev. B 7 , 2336 (1973). 

[18] M.E. Taylor, Ultramicroscopy 42 , 215 (1992). 

[19] M. Schackert, T. Markl, J. Jandke, M. Holzer, S. Ostanin, 
E. K. U. Gross, A. Ernst, W. Wulfhekel, Phys. Rev. Lett. 
114 , 047002 (2015). 

[20] Alan J. Bennett, C. B. Duke and S. D. Silverstein, Phys. 
Rev. 176 , 969 (1968). 

[21] C. Berthod, T. Giamarchi, Phys. Rev. B 84, 155414 

( 2011 ). 

[22] J. Bardeen, Phys. Rev. Lett. 6, 57 (1961). 

[23] M. H. Cohen, L. M. Falicov, and J. C. Phillips, Phys. 
Rev. Lett. 8, 316 (1962). 

[24] Gerald D. Mahan , Many-Particle Physics , Springer, 3rd 
ed. (2000) . 

[25] J. R. Kirtley and D. J. Scalapino, Phys. Rev. Lett. 65, 
798 (1990).' 

[261 Ming-wen Xiao, Zheng-zhong Li, Physica C: Supercon¬ 
ductivity 221 , 136 (1994). 

[27] A. V. Gold, Philos. Mag. 5, 70 (1960). 

[28] L. Zhang, T. Miyamachi, T. Tomanic, R. Dehm, and W. 
Wulfhekel, Rev. Sci. Instr. 82 , 103702 (2011). 

[29] W. McMillan , Phys. Rev. 167 , 331 (1968). 

[30] Brun, C. and Hong, I-Po and Patthey, F. and Sklyadneva, 
I. Yu. and Heid, R. and Echenique, P. M. and Bohnen, K. 
P. and Chulkov, E. V. and Schneider, Wolf-Dieter, Phys. 
Rev. Lett. 102 , 207002, (2009). 

[31] D. Eom, S. Qin, M.-Y. Chou, and C. K. Shih, Phys. Rev. 
Lett. 96 , 027005 (2006). 

[32] I. B. Altfeder, K. A. Matveev, and D. M. Chen, Phys. 
Rev. 78 , 2815 (1997). 










6 


[33] H. H. Weitering, D. R. Heslinga, and T. Hibma, Phys. 
Rev. B 45 , 5991 (1992). 

[34] M. Jalochowski, H. Knoppe, G. Lilienkamp and E. Bauer, 
Phys. Rev. B 46, 4693 (1992). 

[35] I. Yu. Sklyadneva, R. Heid, K.-P. Bohnen, P. M. 
Echenique and E.V. Chulkov, Phys. Rev. B 87, 085440 
(2013) . 

[36] R. Heid, K.-P. Bohnen, I. Yu. Sklyadneva, and E. V. 
Chulkov, Phys. Rev. B 81 , 174527 (2010). 

[37] B. N. Brockhouse, T. Arase, G. Caglioti, K. R. Rao, and 
A. D. B. Woods, Phys. Rev. 128, 1099 (1962). 

[38] Charles Kittel, “Introduction to solid state physics”, Wi¬ 
ley New York (2005). 

[39] T. Nishio, M. Ono, T. Eguchi, H. Sakata, and Y. 
Hasegawa, Appl. Phys. Lett. 88, 113115 (2006). 

[40] T. Nishio, T. An, A. Nomura, K. Miyachi, T. Eguchi, H. 
Sakata, S. Lin, N. Hayashi, N. Nakai, M. Machida, and 
Y. Hasegawa, Phys. Rev. Lett. 101 , 167001 (2008). 

[41] S. Qin, J. Kim, Q. Niu, C-K. Shin, Science 324, 1314 
(2009). 

[42] A. M. Garcia-Garcia, J. D. Urbina, K. Richter, E. A. 
Yuzbashyan, B. L. Altshuler, Phys. Rev. B 83, 014510 
( 2011 ) 

[43] S. Bose, A. M. Garcfa-Garcfa, M. M. Ugeda, J. D. Urbina. 
C. H. Michaelis, I. Brihuega and K. Kern, Nat. Mat. 9, 
550 (2010) 

[44] A.A. Galkin, A.I. D’yachenko, V.M. Svistunov, JETP, 
Vol. 39, No. 6, p. 1115 (1974) 

[45] J. Schmalian, M. Langer, S. Grabowski and K. H. Ben- 
nemann, Computer Physics Communications 93, 141 
(1996). 

[46] J. Rossat-Mignod, L. Regnault, C. Vettier, P. Bourges, 
P. Burlet, J. Bossy, J. Henry, and G. Lapertot, Physica 
C 86, 185 (1991). 

[47] H. A. Mook, M. Yethiraj, G. Aeppli, T. E. Mason, and 
T. Armstrong, Phys. Rev. Lett. 70 , 3490 (1993). 

[48] H. Fong, P. Bourges, Y. Sidis, L. Regnault, A. Ivanov, 
G. Gu, N. Koshizuka, and B. Keimer, Nature 398, 588 
(1999). 

[49] A. D. Christianson, E. A. Goremychkin, R. Osborn, 
S. Rosenkranz, M. D. Lumsden, C. D. Malliakas, I. S. 
Todorov, H. Claus, D. Y. Chung, M. G. Kanatzidis, R. 
I. Bewley, and T. Guidi, Nature 456 , 930 (2008). 

[50] D.S. Inosov, P. B.J.T. Park, D.L. Sun, Y. Sidis, A. 
Schneidewind, K. Hradil, D. Haug, C.T. Lin, B. Keimer, 
and V. Hinkov, Nat. Phys. 6, 178 (2010). 

[51] A. Abanov, A. Chubukov and J. Schmalian, J. Electron 
Spectrosc. Relat. Phenom. 117 , 129 (2001). 

[52] H. Suderow, E. Bascones, A. Izquierdo, F. Guinea, S. 
Vieira, Phys. Rev. B 65 , 100519(R) 2002. 

[53] A. Kamenev, Field theory of non-equilibrium systems, 
Cambridge Univ. Press (2011) 

[54] A. I. Larkin and Y. N. Ovchinnikov, Sov. Phys. JETP 
41 , 960 (1975) 

[55] Note that the extensions of the lead islands are typically 
larger than the 400 x 400 nm 2 STM images of the sur¬ 
faces giving a minimal island size of 0.16 pm 2 . 

[56] Note that also the second peak may already include such 
two-phonon processes. 

[57] Note that we should not broaden the inelastic contribu¬ 
tion as we use the a 2 ftun(ca) from the normal conductor 
measurement that already includes broadening, see Sup¬ 
plementary Material for details. 


Supplemental Material 

Derivation of the tunnel current 


Peturbative approach 


The tunneling current is given by elementary charge 
times the change of the number of electrons ns = 
Sfc a c l < 7 C k,cr in the superconductor 

I = -e-^tr [p(t)n s \ /tr [p(t)} 

= ie tr(p 0 [n s {t),H(t)])/ti[po] 

= ie([n s (t),H(t)\) (5) 


where p(t) = U(t,oo)poU'(t,—oo) is the time-dependent 
density matrix and (... )o = (po ■ ■ ■) = ( ...) is the 

expectation value of the system in thermal equilibrium 
with density matrix p 0 = eA suitable formalism to 
calculate the current © is the Keldysh Green function 
method (we follow the notation of Ref. [53]). The cor¬ 
responding Keldysh action of the Hamiltonian (without 
bias voltage) employed in the main text of the paper is 
given by 


S = So + S t 


( 6 ) 


So = dt - €p )c p ,a{t) 

-' C P,CT 

+ dt ~ e k)Ck,a(t) 

^ C k,a 

+ / dt y ^ - uJ q ,p)a q ^{t) 

Jc 

j I dt \ (t)Cfc' ,cr (t'j^k — k' ,fi (t) 

VVs Jc l 


k,k' 


S, = - 


W^Vi d dt T^pCk,a{t)Cp : g(t) 

7wl c dt Tip ’ q 


Vs\/Vt Jc 
+ h.c. 


[ U^g,/iCfc i(T (t)Cp,cr (t)0g,^ it) 


where as usual we defined the phonon displacement field 
= dq^+aLq ^[? ]. In order to derive the I—U char¬ 
acteristic of the system we have to apply a finite voltage 
eU = ps ~ Hr, which can be done easily by the substitu¬ 
tion Cp, CT (f) -A e -l/iTi c Pi<T (f) and c fcj(T (f) -4 e -lMst c fci<J (f). 
The dispersion energies of the tip — pt 

and superconductor -A = ef — pr are now mea¬ 
sured relative to their chemical potentials and this leads 
to a time dependence of the tunneling matrix elements 
T e -A T e e leUt ,T 1 -A in the tunneling part S t of 

the action. 

We build up our perturbation theory by rewrit¬ 
ing (...) = J D[c, c]D[d, a].. ,e lS = (e lSt ...) 0 

with the unperturbed expectation value (■ ■ ■ )o = 
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f D[c, c]D[a , a]... e lS °. The corresponding unperturbed 
propagators are then given in the i?, A, K basis (also 
known as Larkin-Ovchinnikov Representation) 



Gk/ P {t, t') 


fe(M') 




f C k/p,a 

\ C k/p,a 



( C k/p,a 

\ C k/p,a 



G^ /v (t,t') 


~l{ 


0 


k/p' 

G£/p(t,t') y 

feY )o 

fej )0 


fefe') D*Jt,t')\ 

Vfe(M') 0 ) 


)o 


( 7 ) 


FIG. 6: Normal and anomalous self-energy due to electron- 
phonon interaction that appear in the Eliashberg-theory. 


the Keldysh propagators have the simple structure 
G % p M = [1 - 2n F (u)\ - fe(W] , 

s -V-' 

-2mA k/p {w) 

= [1 + 2 n B (co)\ - fe»] , (9) 

-2iri [a 9iM (uj)— w)] 


where the retarded propagators are given in energy rep¬ 
resentation f(w) = f^o dtf(t)e lult as 


Gfc(w)-•— 

feW = _ 

Gp (w) = - ♦ - 

CM = wv 


_ -feW _ 

[Zf( W ) W ] a -[ff] 2 -$»> 


l 

w - + iO 

(w + *0) 2 - fe. 


( 8 ) 


where n F {<jS) is the Fermi and ns(w) the Bose function 
with temperature T and we defined the spectral weights 
Ak/p{u}), A qt p(ui) of the electron and phonon systems. In 
our case A q fJj (uj) = S(oj — ui q p). For completeness, we 
also give the explicit expressions for the greater/lesser 
Green functions 



( 10 ) 

A qA~ U ) 


Following Eq. © it is easy to determine the explicit ex¬ 
pression for the current 


For the superconductor, we use the known framework of 
Eliashberg theory. Therefore, Zg = 1 - ^G)-^- k (-^) 
is the renormalization function of the lead superconduc¬ 
tor S and = Efc(w ± i 0) is the phonon-induced 

normal self-energy, see Fig. [G] We neglected the correc¬ 
tion of the pure dispersion due to the coupling to 
the phonons, because it will basically just give an unim¬ 
portant shift of the chemical potential and can be as¬ 
sumed to be incorporated in the electronic dispersion al¬ 
ready from the beginning. The anomalous self-energy 
(w) = <hfc(w + i0) is depicted in Fig.jbjand we also gave 
the expression for the anomalous propagator -F^(w), even 
tough we will only need the normal particle propagator 
(since in the NIS-junction there is no Josephson effect). 
We neglect the renormalization of the phonon spectrum 
due to the interaction with the electrons, which could be 
incorporated easily by a phonon self-energy that would 
just lead to a small broadening and modification of the 
phonon spectral weight. The effect of the Coulomb in¬ 
teraction between the fermions is as usual incorporated 
using a Coulomb pseudopotential p *. Since we consider 
the sub-systems S and T to be in thermal equilibrium, 


I = 


ie fe( 

k,p 

ie fe( 

k,p 

-«<E 


feW y 
few; 1 

( 0 
few 

-fewr^ 

(cpAt)\ \ 

few; } 

fewV i 

( 0 

-fewn 

( C p,cr(t)\ iSt \ 

few; 1 

few 

o ) 

few; e ' 

AV.crWV 


-few n 

(tfe 

few; 

fe P w 

o ) 


(ii) 


where we defined the total tunneling matrix ele¬ 
ment fe(i) = e ieUt [T| p + E( wi t h 
phonon field A on the upper/lower Keldysh contour) 
and in the end expanded in leading order of T^.p. Also, 
we defined the creation operators to be defined on the 
lower Keldysh contour (index —) and the annihilation 
operators on the upper Keldysh contour (index +), since 
the electron first has to leave one side before it can tun¬ 
nel through the barrier to the other side. This time or¬ 
dering can be conveniently expressed in Keldysh theory 
and we also assumed the phonons (field A) to be ex¬ 
cited/absorbed on the same contour as the electrons on 
the superconductor S. The tunneling action can be writ- 

















ten in a similar way as 


near the Fermi surface, we can then rewrite the elastic 
current (14) to be 





fCp,rr(t)\ 
\Ck,cr{i) J 



P£p(t)] 

0 



( 12 ) 


Elastic current 


Performing the contractions in Eq. @ we find the 
elastic current 


r((/)=2e i dT vk ^ ,Ti ' 

k,p o:=i 


I ^ gieUr 


E 


[cr (-t)G+’“(t) - G+-“(—t)G“-(t)] 

oo 

1 


= —4e / dr 


VsVt 


EI t *, ? 


12 ^ieUr 


(13) 


k,p 


[im Gfc (—r)Im G^ (r) - ImG^(-r)ImG<(r)] 


where we used the definition of the greater/lesser 
G > / < and time-ordered/anti-time-ordered Green func¬ 
tion G r ^, see Section | of the Supplementary Material. 
In the end, we used the known identities that relate G T ^ 
to the greater, lesser, retarded and advanced propaga¬ 
tors. In Fig. [7] we show the corresponding Feynman di¬ 
agram for the elastic tunneling current for the leading 
order in the tunneling element. After transforming to 
Fourier space and inserting the explicit expressions of the 
propagators in thermal equilibrium, we find 


OO 

nu>=i*e J 

-OO k ’P 


(14) 


[n F (uj) — n F (w + eU)\A k (w)Ap(w + eU), 


which is the usual expression for the elastic current in 
the Landauer-Biittinger transport theory assuming per¬ 
fect quasiparticles A k / p (oj) = 5 (lj — e^J). 

The elastic conductance G e (U ) = ^ will now be cal¬ 
culated using the usual approximation T% = t e = const, 
for small voltages U <C Ep ~ 1 eV, which is a reasonable 
assumption for an STM [23] Assuming the DOS of the 
tip system i/ T (o;) = 1/Vr Y P ^p( w ) ~ u t t° be constant 


I e (U) = 4t rz/£e 



dui [n F (u>) 


n F (u + eU)]u s {u }), 



FIG. 7: Feynman diagrams for the elastic (left) and inelastic 
(right) tunneling current in leading order T e ,T\ 


where we defined as usual the DOS of the superconduc¬ 
tor as v{u)) = 1 /Vs Yk ^fc( w )- As it is we ll known, the 
elastic differential conductance is then given by 




du)n' F (u) + eU)us{ui) 

J — OO 

/ OO 

duj n' F (co + eU)v , s{ijj) (16) 

-OO 


where we defined as the normalized DOS v(uj) = v(oj)/v^ 
with i/g as the DOS and op = 47r^^e 2 |t e | 2 as the elastic 
conductance in the normal state. For small temperatures 
(T -C E F in the normal conductor or T « A in the 
superconductor with gap A), such that n' F (e) ~ —6(e), 
the elastic conductance simplifies to 


dT e 

= 4nv F e 2 \t e \ 2 v(—eU) = a 0 v(—eU ), 


(17) 


and is then proportional to the normalized DOS v(co) 
of the superconductor. The corresponding expression 
von d 2 I e /dU 2 can be computed easily from the expres¬ 


sion (16). 


Inelastic current 

The inelastic current can similarly be expressed by per¬ 
forming the contractions of @ containing the phonon- 
fields and expressing the occurring propagators in terms 
of retarded, advanced, greater and lesser propagators 
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00 

\U) = 2ie f dr ' £ e ieUT £ a [G“--(-r)G+’“(r)D-“(r) - G+’“(-t)G“>-(t)D“;+(t) 

5 T a=± L 

9,A 4 


— OO 

oo 


= 4e / dr 


V * I 2 V T 

<0 /c,p 


\ ' 17^1 


I 2 gief/r 


Im (-r)Im G< (r)Im D^(t) + Im (-r)Im G~(r)Im (r) 
-ImGf(-r)ImG<(r)Im£i< /i (r) (18) 


After going to Fourier space and inserting the corresponding electron and phonon propagators defined in Sec. [] we 
can finally rewrite the inelastic current as 


P(U) = —47re J dujiduj2 


' v§v T 


EK 


p , q , p Q , M 


( 19 ) 


A g ,^(cJi)A fc (w2)Ap(a;2 - wi + eU) yn F (oo2 - u>i + eG)ras(wi)[l - n F (w 2 )] - nj?(w 2 )[l + n B {ui)\ [l - n F (w2 - wi + eG)] J 
-t- Ag ! ^ i (o;i)Afe(aJ2)Ap(a’2 -t un + eU) I tlf(oj2 + cai + eG) [l + tib(m 1)] [l — ^^(^2)] — tif (<^2)^.23(021) (l — w F (w2 + un + eG)] 



FIG. 8: Inelastic tunneling processes for F > 0 and V < 0 
via the emission/absorption of phonons. 


The first/third term describes the tunneling of an elec¬ 
tron from the tip to the superconductor via the absorp¬ 
tion/excitation of a phonon and the second/fourth term 
the tunneling from the superconductor to the tip via a 
phonon excitation/absorption, see also Fig. [8] As in the 
elastic case, we apply the following simplifications: A 
constant inelastic vertex = t l and a constant 

DOS of the tip. Let us define the weighted DOS of the 
phonons in the superconductor as 

a Ftmi(w) = y~ l a q, M | A q ^(u}) 

9 

= y~ l a g,A»l — Wg,/i) > (20) 

' q,F 

which is very similar to the Eliashberg function besides 
a different momentum average. For our case of very low 
temperature ksT <C ojd only the processes that excite a 
phonon lead relevant inelastic contributions to the tun¬ 
neling current since the number of thermal low-energy 
phonons a 2 Ft un (w) • «s(w) ~ 0 in the system is negligi¬ 


ble. We then find 


I 1 = 4irei' F Vg t* | / doJidu} 2 ct -Ft un (wi)^ s (w 2 ) (21) 


n F (u) 2 )n F (uJi — w 2 — eU) - {w 2 , U —► —w 2 , —£/} 

For the differential conductance we then find 

2 „ 


dP_ 

dU 


= -CTq 



k 

2 f 



J [ 


duj\ dui2 ck 2 Etu„ (wi) Fs (w 2 ) (22) 

n_F(w 2 )n' F (wi -u) 2 - eU) + {w 2 , [/ ->■ — w 2 , —J7} 

« 2 ^ T u„(- + eU)vs(w)n F {u) 

+ a 2 Ft un {u-eU)v s {-u)n F (u) (23) 


where we defined the thermal broadened weighted DOS 
of the phonons as the convolution (in the limit of zero 
temperature it obviously holds a 2 F^°(uj) = F tnn (uj)) 

/ OO 

dy a 2 F tun (y)n' F {y - x) (24) 

-OO 

For particle-hole symmetric electronic systems f's(w) = 
—vs(—io) and in the limit that k F T is much smaller than 
the characteristic phonon frequencies (meaning F^ ln (w ± 
eU)n F (uj ) ~ 0(=F U)) , we can simplify this expression to 


dF 

dU 


= a 0 


J dcua 2 F^ un (uj + e\U\)i>s{uj)nF(uj) ( 25 ) 
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FIG. 9: Elastic (blue dashed), inelastic (green dashed) con¬ 
tribution to the total (black solid) second derivative of the 
current d 2 I/dU 2 in the superconducting state for different 
peak width 70 . The additional inelastic contribution will lead 
to a peak-dip feature with similar positive and negative am¬ 
plitude in the second derivative, whereas the purely elastic 
contribution only pronounces the dip at the shifted phonon 
frequency eU as A + ujq strongly. 


Elastic and inelastic STM tunneling for single mode 


In order to get a qualitative understanding of the 
inelastic tunneling contribution in the superconducting 
state, let us analyze a simple toy model. The toy model 
consists of a single phonon mode with a 2 F tnn (u>) ~ 
a 2 F(uj) = A 0 f(oj) ( ^_^ 2 +7 2 with characteristic phonon 
energy wo = 5 meV and half-width 70 . The function 
/(w) = u 2 +(i meV) 2 ensures the proper low frequency be¬ 
havior of acoustic phonons and rapidly approaches unity 
for larger frequencies. The value A 0 is choosen such 
that the dimensionless electron-phonon coupling constant 
A = 2 f^° du a 2 F(u)/uj = 1.5. We use fi* = 0.1 for the 
pseudopotential, such that solving the Eliashberg equa¬ 
tions yields a gap value A ~ 1 meV. In Fig. [9] we see 
the resulting second derivative of the tunneling current, 
which is more sensitive to the fine-structure than the con¬ 
ductance, for the above mentioned elastic and inelastic 
tunneling for different peak width 70 in the supercon¬ 
ducting state. As was seen in the experiments we use 
the ratio | i 1 /| 2 as 0.12/ meV duj a 2 F(u). The inelas¬ 
tic contribution Eq. ( |22[ ) lias its major contribution for 
frequencies a bit below the energy of the phonon peaks 
shifted by the gap A. Since inelastic tunneling adds addi¬ 
tional channels to the conductance, its contribution will 


have the opposite sign of the elastic contribution in (16) 


and can give pronounced positive peaks in the second 
derivative of the tunneling current followed by a nega¬ 
tive peak of same amplitude . These symmetric peak-dip 
features around the zero axis in the second derivative of 


the tunneling current are characteristic for the joint elas¬ 
tic and inelastic STM. In the present example, we were 
not able to find an Eliashberg function that would yield 
such a tunneling spectra from the purely elastic tunnel¬ 
ing contribution Eq. (16). Even for very sharp Eliashberg 
spectra where the second derivative of the tunneling cur¬ 
rent can be clearly positive for some voltages, see upper 
picture in Fig. [9j the following dip will always be much 
more pronounced if one only considers the elastic tunnel¬ 
ing contribution. 


Modelling of experimental broadening 


The experimental resolution is limited due to the used 
lock-in technique. As was shown in Ref.ITZi for a modula¬ 
tion voltage Vjnod (the modulation voltages stated in this 
work are root mean square values U mo a = t/mod.max/v^) 
the experimental curve for the second derivative mea¬ 
surements is the actual current convoluted with a Gaus¬ 
sian function r(w) of half-width FWHM« 1.7ekjnod- 
This corresponds in our case to a standard deviation 
a = FWHM/2.4 = 0.31meV. For the elastic part this 
means that we get 


d 2je,exp 

dU 2 


/ OO 

dEdujT(eU — E)n' F (uj + E)0s(w) 

-OO 

(26) 


which can be easily computed using the fermionic DOS 
obtained from solving the Eliashberg equations. 

Let us now consider the experimental broadening for 
inelastic tunneling and restrict us to the case of positive 


bias voltages. Following Eq. (25), the second derivative 


of the inelastic tunneling current is for positive U > 0 
given by 


d 2 P 

W 2 


= eu 0 


/ ’ OO 

dhja 2 F^{eU + w)i> s (w)n F (w) 

-OO 

(27) 


To get the experimental data (both in the normal and 
superconducting state), we have to broaden the function 
by the convolution 


d 2 P’ exp 
dU 2 


= ecr 0 


= ea 0 


dET(eU-E) 


(28) 


dus a 2 F 2 un '(E + u)v s (u)n F {u) 

OO 

J du a 2 F™ pl (eU + u)D s (uj)n F (u) 

with the thermally and modulation voltage broadened 
spectral function 

/ OO 

dyT{x - y)a 2 F^ un {y) (29) 

-OO 
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In the normal state r's(w) ~ 1 this simplifies in the limit following tunneling processes 
of small temperatures T <C uiD,Ep (such that ?zf(w) ~ 

6{— w)) to 


d 2 J‘’ exp 


dU 2 


ecr 0 


ecr o 


,i - roo 

- / dwa 2 JCn P, (e^ + w)riFH 

^ J —OO 

^E 2 F^(et/) (30) 


As P ~ 0, we can extract the a 2 F^™ od (eU) func¬ 
tion from the normal state measurements and can then 
use it to calculate the inelastic current in the supercon¬ 
ducting state. 


8H[ n) = 


Vv t v< 


^ E 


T 1 


k,P,q l, 
cr,Hl 


• • • ^q n ,fJ-n^k CT^P i^t Qiifl 1 




(31) 


Multiple phonon processes 

If we consider tunneling processes with a higher num¬ 
ber of excited phonons we can formally write down the 


In the zero temperature limit it is then straightforward 
to generalize the result (21) to the n-phonon process (de¬ 
manding energy conservation and Fermi statistic for the 
leads) 


/Un) = 47rel ,o 




J duj i ... dcu n duj n -\-iQz Ft un (uq)... Oi Ft un (w n )F5(c^ n -)_i) 


(32) 


Lo n+ i)6{uj n+ i — iu n — ... — u )i + eU) — {w n+ i, U —u> n+ i, —U} 


For the conductance, we then find for particle-hole symmetric systems 


dl' 1 '^ 

dU 


= (TO 


= <TO 


t 0{n) 


sign (U) / dioduj!... du n -iot 2 F tun (e \U\ - u - U}i)a 2 F tun (u)i - w 2 )... a 2 F tun (uj n - 2 - w„_i)i>s(w) 


t T(n) 


t e 


sign( U) / dLua 2n F£ an (e\U\-u)i>s(uj) 


)Q 


where we defined the convolution 

_ 2 n rin 


a 


/*oo 

-^tun(^) = / did\ . . . dcd n — i OL CJi)oi -ftun(^l ^2) • • • ^ -^tun(^n—2 ^n— l) 

JO 


(33) 


(34) 


Important relations of Non-Equilibrium propagators 


We can now perform the Keldysh rotation to the classical 
and quantum fields in the bosonic case: 


Following Ref. [53] we define for both the fermionic and 
the bosonic fields <j>(t) the greater, lesser, time-ordered 
and anti-time-ordered Green’s functions as 


G T (t,t') = G ++ {t,t') = -i{cj> + {t)4> + (t')) 
G f it , if) = G (t, t') = -i{ (j)~ (t)<j>~ {t')) 


+ 0 (*)] 

W)E[^(t)-r(i)] 


(36) 


and similar for the conjugate fields <j>{t). However, for 
the fermionic fields we use the Ovchinnikov-Larkin con- 
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vention [M] 


frit) = ^[</> + W + <TW] 
frit) = ^|> + (*) - <P~(t)\ 
frit) = ^[4> + (t) - 4>~{t)\ 
frit) = -^\4> + {t) + 4>~{t)\ 


(37) 


For the fermionic and bosonic cases, the retarded, ad¬ 
vanced and Keldysh propagators are then defined as in 


Eq.fTf. The relations between the different Green’s func¬ 
tions (>, <,T,T, R, A, K ) can be summarized by 


0 = G r + G f - G> - G< 

G K = G> + G< 

G r = \ [G t - G' r + G > - G<] 
G a = \ [G t -G f -G> + G<] 


